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00 , The quantum integrabiUty is established for the one-dimensional supersymmetric U model with 

0^ ■ boundary terms by means of the quantum inverse scattering method. The boundary supersymmetric 

U chain is solved by using the coordinate space Bethe ansatz technique and Bethe ansatz equations 
are derived. This provides us with a basis for computing the finite size corrections to the low lying 
M ' energies in the system. 
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In recent years, there has been a considerable interest in exactly solvable lattice models vi/ith boundary fields and/or 
interactions |^-|^. One class of such models are one-dimensional boundary strongly correlated electron systems, which 



I is of great importance due to their promising role in theoretical condensed-matter physics and possibly in high-Tc 



superconductivity Q. Boundary conditions for such systems, which are compatible with integrability in the bulk, 
^ . are constructed from solutions of the (graded) reflection equations (called boundary K-matrices) | i[| . Work in this 
direction has been done for the Hubbard-like models (|-|^ and for the supersymmetric t-J model p^pT |. 

In this Letter, we study integrable open-boundary conditions for the supersymmetric U model of strongly correlated 
electrons introduced in |l^ , |l3| and extensively investigated in [^ - [l6| . We will present a boundary supersymmetric 
''O ' U model and show that it can be derived from the quantum inverse scattering method by modifying and generalizing 
C ' Sklyanin's arguments, thus establishing the quantum integrability of the boundary model. In doing so, we encounter 
the following complication: the Sklyanin's definition for boundary Hamiltonian can not apply since the supertrace 
of the boundary K-matrices of zero spectral parameter is equal to zero for the current case. This is related to 
the fact that the supersymmetric U model has been constructed from the R-matrix associated with the typical 4- 
' dimensional irreducible representation of gl{2\l). Nevertheless we manage to solve this complication by introducing a 
new definition for Hamiltonian. We then solve the boundary supersymmetric U model by the coordinate space Bethe 
\^ ' ansatz approach and derive the Bethe Ansatz equations. 

. Let c| ^ and Cj.a- denote fermionic creation and annihilation operators with spin cr at site j, which satisfy the anti- 

commutation relations given by {c| ^, Cj,r} = SijSa-r, where i, j = 1, 2, • ■ • , L and cr, r =|, j. We consider the following 
Hamiltonian with boundary terms 



L-1 



td' H=J2Hfj+,+Bu + Bru (1) 



' O ' where Hf-,^ is the local Hamiltonian of the supersymmetric U model introduced in pj 



cr 

+ ^(^31^31 + "j+iT%+u) + *p(4TsVj+UCi-l-iT + h-c.) + (rij + n^+i), (2) 



and Bit, B^t are boundary terms 



2{U + 2) (2 \ u 2(;7-H2) (2 
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In the above equations, nj„ is the number density operator njcr = Cj^cja, nj = riji and tp = ^ = ±(1 — exp(— 77)); 
^± are some parameters describing boundary effects. 

Some remarks are order. As is seen from (|^), Bn (Brt) is an inhoniogenous combination of two terms contributing 
to the left (right) boundary conditions. The physical meaning of these terms in the context of strongly correlated 
electrons are the following. The first term is nothing but a boundary on-site Coulomb interaction and the second 
term is a boundary chemical potential. 

We will establish the quantum integrability of the boundary supersymmetric U model (|l|) by showing that it can 
be derived from the quantum inverse scattering method. Let us recall that the Hamiltonian of the supersymmetric 
U model with the periodic boundary conditions commutes with the transfer matrix, which is the supertrace of the 
monodromy matrix T{u), 

T{u)^ Rol{u)---Roi{u). (4) 

The explicit form of the quantum R- matrix Roj{u) is given in p^ . Here u is the spectral parameter, and the subscript 
denotes the auxiliary superspace V — C^'^. It should be noted that the supertrace is carried out for the auxiliary 
superspace V. The elements of the supermatrix T{u) arc the generators of two associative superalgebra A defined by 
the relations 

i?i2(ui - U2) T (wi) T (W2) =T (W2) T (Mi)i?i2(wi - U2), (5) 
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where X= X ^1, X= 1® X for any supermatrix X G End{V). For later use, we list some useful properties enjoyed 
by the R-matrix: (i) Unitarity: Ri2{u)R2i{—u) = 1 and (ii) Crossing-unitarity: i?^2^ (— u + 2)i?2i^ (it) = plu) with 
p{u) being a scalar function, p{u) = v?'{2 — /[{2 + 2a — u)'^{2a + u)"^]. Throughout this letter, a = ^. 

In order to describe integrable electronic models on open chains, we introduce an associative superalgebras and 
7^ defined by the R-matrix R{ui — U2) and the relations 

1 22 1 

i?12(ui - U2) T- (Ml)i?2l(wi + U2) T- {U2) =T- {u2)Rl2{ui + U2) T- (ui)i?2i(wi ~ U2) , (6) 



Rf/''H~Ui+U2)T+ {ui)Ri2{~Ui-U2+2)T+ {U2) 

=r+ ' {u2)R2i{~ui -U2 + 2)T+ iui)Rll''''^ i-ui + U2) (7) 

respectively. Here the supertransposition st^ {fi — 1,2) is only carried out in the /i-th factor superspace oiV(!S)V, 
whereas istf^ denotes the inverse operation of sf^. By modifying Sklyanin's arguments |jl[, one may show that the 
quantities t(m) given by t{u) = str{T+{u)T-{u)) constitute a commutative family, i.e., [t{ui),t(u2)] — 0. 
One can obtain a class of realizations of the superalgebras T+ and 71 by choosing T± {u) to be the form 

T^{u) = T^{u)f^{u)Tl\~u), Tfiu) = Tfiu)ff{u) {T^\~u)f (8) 

with 

T^{u) = Rqm{u) ■ ■ ■ Roi{u), T+{u) = RoL^u) ■ ■ ■ Ro,M+i{u), 'f±{u) = K±{u), (9) 

where K±(u), called boundary K-matrices, are representations of T± in Grassmann algebra. 

We now solve (^) and (0) for (u) and if _ (u) . For simplicity, let us restrict ourselves to the diagonal case. Then, 
one may check that the matrix K-(u) given by 

= ^_(2^_g_) di^g {A^{u),B^{u),B^{u),C^{u)), (10) 

where A_{u) = (^_ + u) (2 - ^_ - u) , S_ (u) = (^_ - u) (2 - ^_ - u) and C_ (u) = (C- - u) (2 - ^_ + it) , satisfies (|) . 
The matrix K^{u) can be obtained from the isomorphism of the superalgebras 71 and 71^. Indeed, given a solution 
T_ of (|), then T+{u) defined by 

Tf{u)^T^{~u + l) (11) 

is a solution of (|^). The proof follows from some algebraic computations upon substituting (11) into (^ and making 
use of the properties of the R-matrix . Therefore, one may choose the boundary matrix K^(u) as 
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i^+H =diag {A+{u),B+{u),B+{u),C+{u)) (12) 

with A+{u) = (2 - 2a - C+ - u){2a + C+ + B+{u) = {-2a - ^+ + u){2a + ^+ + u) and C+{u) = (-2a - C+ + 
u){2 + 2a + ^+-u). 

Now it can be shown that Hamiltonian (|l|) is related to the transfer matrix t(u) (up to an unimportant additive 
constant) 

^^_2_(U+^ 
U 



4:{V + 2W) ^ 2 ' ' 2(y + 2T4^) 

+2stro (k'+ (0)i?fo ) + stro (k+ (0) (i?fo) 



stro ( K+ (O)Glo 



where 



(13) 



V = stroX;(0), = stro iyK+ (O)i/fo y > 

H^^ = /^,,i?^,,(0), G,,, = P^,JK,(fi). (14) 



Here j denotes the graded permutation operator acting on the i-th and j-th quantum spaces. ( p_3D imphes that the 
boundary supersymmetric U model admits an infinite number of conserved currents which are in involution with each 
other, thus assuring its integrability. It should be emphasized that Hamiltonian (|l|) appears as the second derivative 
of the transfer matrix t(u) with respect to the spectral parameter u at u = 0. This is due to the fact that the 
supertrace of K-^-{0) equals to zero. As we mentioned before, the reason for the zero supertrace of -^^4.(0) is related to 
the fact that the quantum space is the 4-dimensional typical irreducible representation of gl{2\l). A similar situation 
also occurs in the Hubbard- like models 

Having established the quantum integrability of the boundary supersymmetric U model, we now solve it by using 
the coordinate space Bethe ansatz method. Let us assume that the eigenfunction of Hamiltonian ^ has the form 

N 

'^ai,- -,aN{xi,- ■ ■ ,xn) = ^ ep^<TQi ,...,o-q„ (fcpQi , ' ' ■ , kpQN) exp{i^ kp^ X j) , (15) 

p j=i 

where the summation is taken over all permutations and negations of fci , • ■ • , , and Q is the permutation of the N 
particles such that 1 < xqi < ■ • • < xqn < L. The symbol ep is a sign factor ±1 and changes its sign under each 
'mutation'. Substituting the wavefunction into the eigenvalue equation H\'i>) = E\'i>), one gets 

^ - ■ ■ ,(Tj ,cri , - - ■ (■ ■ ' 1 kj ^ ki^ • • ■) — Sij {ki ^ kj^ A.. . . fj^ fjj . • • ^ ki^ kj J • • • ) , 
A^^,...{-kj, • ■ •) = s'^ikj;pia,)A^^,...{kj, ■ ■ •), 

A..,,, (• ■ • , -fcj) = s^ikj;pLa.)A....,,, (■ • ■ , fcj) (16) 
with Sij{ki, kj) being the two particle scattering matrix and s^, the boundary scattering matrices, 

e{ki) - e{k,) + icv.j 



6{ki) - 6{kj) + ic 
1 - Pia, e.'^' 



where pi„ = pi = — 1 + ^^'^^'^ 2-; ' ^'^'^ = Pl = ^1 + 2-;+ ^^"^ ^ e'' — 1; "P^ is a spin permutation 

operator, and the charge rapidities 0{kj) are related to the single-particle quasi-momenta kj by 9{k) = itan(-|) [Q. 
Then, the diagonalization of Hamiltonian reduces to solving the following matrix eigenvalue equation 
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Tjt = t, j = l,---,N, (18) 
where t denotes an eigenvector on the space of the spin variables and Tj takes the form 

T, = S-{k,)s''{-kj;pi„^)Rj{k,)R+ikj)s''{k,;pLa,)S+{k,) (19) 

with 

Sj'ikj) = Sj^Nikj,kN) ■ ■ ■ Sj.j+i{kj,kj+i), 

(%) ~ '^i,i-i(^i> ' ' ' ^j-ii^j^ ^i)' 
(%) = Sij{ki, —kj) ■ ■ ■ 5j_ij(fcj_i, —kj), 
R^{kj) = Sj+i.j{kj+i,-kj) ■ ■ ■ SN.j{kN, -kj). (20) 

This problem may be solved using the algebraic Bethe Ansatz method. The Bethe ansatz equations are 

M 



Ag - 6j + if Ag + 6j + _ -pj Ag - A/3 + jc Xa + \0 + ic 

where 6j = 6{kj) and C{k;p) = (1 — pe^''^)/(l — pe**^). The energy eigenvalue E of the model is given by i? = 
— 2^^.^-^ cos kj, where we have dropped an additive constant. 

In conclusion, we have studied integrable open-boundary conditions for the supersymmetric U model. Its quantum 
intcgrability follows from the fact that the Hamiltonian of the model on the open chain may be embbeded into a one- 
parameter family of commuting transfer matrices. Moreover, the Bethe Ansatz equations are derived with the use of 
the coordinate space Bethe ansatz approach. This provides us with a basis for computing the finite size corrections 
to the low lying energies in the system, which in turn allow us to use the boundary conformal field theory technique 
to study the critical properties of the boundary. The details will be treated in a separate publication. 
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